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The confined variational method is used to generate a basis of correlated gaussians to describe 
the interaction region wave function for positron scattering from the H2 molecule. The scattering 
length was » —2.7 ao while the zero energy Z e a of 15.7 is compatible with experimental values. The 
variation of the scattering length and Z e « with inter-nuclear distance was surprisingly rapid due to 
virtual state formation at R « 3.4 ao. 
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The lack of spherical symmetry makes the calculation 
of electron or positron scattering from molecules an es- 
pecially intractable computational problem. The non- 
spherical potential couples different partial waves result- 
ing in an enormous escalation in the size of the calcu- 
lation when compared with atomic targets. One conse- 
quence of this is that it is difficult to identify a defini- 
tive calculation of low energy electron/positron scatter- 
ing from the simplest of molecules, i.e. H2, even under 
the simplifications of the fixed nucleus approximation. 

A new approach to compute the wave function for 
electron/positron scattering from small molecules is de- 
veloped. It utilizes existing computational technologies 
from few-body physics that had been used to describe 
the low energy scattering of simple and composite pro- 
jectiles from atoms [H, [H y] . The method is applied to the 
calculation of positron scattering from the H2 molecule. 
The cross section for positron annihilation at thermal 
energies was found to be compatible with experimental 
values [1, H, @ • This is a significant achievement since 
the annihilation cross section presents a stringent test 
to the accuracy of the scattering wave function 0] and 
its successful prediction solves a previously intractable 
problem. Our calculations also show the existence of an 
unexpected virtual state at a H2 inter-nuclear distance 
of i?« 3.4 o . 

There have been a number of calculations of low en- 
ergy e + -H 2 scattering and annihilation 0, H, OH [H| • 
At present, all previous calculations significantly under- 
estimate the low energy annihilation cross section. The 
most sophisticated calculations are the Kohn variational 
calculations performed by Armour and colleagues at the 
University of Nottingham (UN) [1 [O, [ll| . Their most 
recent calculations significantly underestimate the anni- 
hilation cross section at thermal energies. 

We apply a variant of the confined variational method 
(CVM) [J [4] to describe low energy positron-H.2 scatter- 
ing. In the CVM, an artificial confining potential is added 
to the scattering Hamiltonian thus converting the system 
into a bound system. This provides a framework that 
permits the wave function in the interaction region to be 



obtained with bound state techniques. Of crucial impor- 
tance to this exercise is the use of the stochastic varia- 
tional method (SVM) [H, [H, [3 to describe the inter- 
action region wave function. The SVM and variants [l6| 
constitute a powerful tool for studying few body systems. 
The SVM uses a wave function that is a linear combina- 
tion of explicitly correlated gaussians (ECGs) which have 
easy to evaluate Hamiltonian matrix elements [3, [l7| - 
Therefore it is feasible to optimize the non-linear param- 
eters of the basis stochastically. Application to molecular 
systems is easy and ECGs have been recently used to de- 
scribe the wave functions of a number of small molecules 
to high accuracy [l8| . The close to zero energy scattering 
parameters were extracted from the interaction region by 
a stabilization technique [l[ and a technique based on the 
energy [l9j . 

The calculation of the interaction region wave function 
proceeded in a manner that was very similar to previous 
ECG based calculations on collision systems [H,13|. The 
Hamiltonian for e + H2 scattering was 
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The positron coordinate is ro while ri and r 2 are the 
electron coordinates. The vector R/2 is the displacement 
of the two protons from the mid-point of the molecular 
axis. The confining potential Wcp(t) has the form 



Wcp(r) = G(r - R a ) 2 Q(r - R Q ) 



(2) 



where 0(r — Rq) is a Heaviside function and G is a small 
positive number. 

The first stage of the diagonalization of Eq. {T]) was to 
use the SVM to generate an interaction region basis of 
energy optimized ECGs. The ECGs were a generalization 
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TABLE I: The convergence of the various properties of the e + -H2 system for the E 3 symmetry at R = 1.4 ao as a function of 
the number of ECGs, N. The first number in the N column is the dimension of the inner region basis while the second entry 
is the dimension of the outer region basis. The energy of lowest energy state in the confining potential is given by the -Bjv 
column. The wave number, k (in a^ 1 ) is that of the lowest energy pseudo-state when the entire basis was diagonalized without 
the confining potential. The scattering length, A scat (in ao) and Z e g were derived from the wave function projections parallel 
(||) and perpendicular (_L) to the inter- nuclear axis, and from the system energy using the soft-box radius (SB). 
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of those used previously in purely atomic calculations 
[l6|. Their functional form was 

4> k = P exp ^-^y^&fc,ij|r t - Sfc^| 2 ^ 

x ex P ^ a Mjl r i- r il 2 ■ ( 3 ) 

\ i=0 j=i+l / 

The vector Sfe,i displaces the center of the ECG for the 
ith particle to a point on the inter-nuclear axis. This en- 
sures the 3-particle wave function is of £ symmetry. The 
values of dk,ij, bk,ij and Sk,i are adjusted during the op- 
timization process. The operator P is used to enforce E g 
symmetry. Each ECG has a total of nine stochastically 
adjustable parameters. 

Table Q] lists the energy of the confined e + -H 2 sys- 
tem for a succession of basis sets. These energies were 
generated with the confining potential parameters G = 
1.55 x 10 -4 and Rq = 18.0 ao- The inter-nuclear separa- 
tion was set to 1.40 ao which is very close to the position 
of the minima in the H2 potential curve. 

Extracting scattering information requires embedding 
the interaction region wave function into a formalism for 
e -H2 scattering. However, one of our major aims is 
to demonstrate that ECG technologies make it easy to 
get a good description of the e + -H 2 collision dynamics. 
Accordingly, attention is focussed on the very low energy 
region where the outgoing wave is essentially spherical. 

There are two advantages to restricting the current cal- 
culation to very low energy. First, the most reliable ex- 
perimental information comes from traditional positron 
annihilation experiments using thermal positrons that 
yield annihilation cross sections at very low energies [2(| ■ 
Second, the collision can be treated as s-wave scattering 
and thus the molecular aspects of the asymptotic wave 
function can be neglected with minimal error. 



Positron annihilation cross sections are reported as 
Z e g, which is interpreted as the number of electrons avail- 
able for annihilation. The annihilation cross section and 
Z e g are related by the identity 

Z e ${k) = , (4) 

where c is the speed of light. In the first Born approxima- 
tion, the number of electrons available for annihilation is 
equal to the number of electrons in the molecule. 

The scattering length and near zero energy Z e f{ were 
extracted from the wave function using a stabilization 
technique [l|. Initially, the energy optimized interaction 
region ECG basis is supplemented by a set of basis func- 
tions to describe the long range part of the e + H2 wave 
function. The functions were 

^uout = V H2 (ri,r 2 )V^(r ) 
ipi(r ) = P exp (~^ ai rl)^ . (5) 

The target wave function, V H2 (i"i,r 2 ) is represented by 
a linear combination of ECGs. A basis of dimension of 
120 gave an energy of —1.17447554 a.u.. The H 2 energy 
at an inter-nuclear separation of 1.40 ao is —1.17447571 
a.u. [Hj|. Our wave function recovers 99.996% of the 
correlation energy of 0.04084 Hartree [2l|. The ipi{ r o) 
are designed to describe the positron at asymptotic dis- 
tances. The afc were an even tempered set given by the 
identity ay = ai/T^ 1 with ai = 18.59 and T = 1.435. 
A total of 36 long range basis functions were added to 
interaction region basis. 

The Hamiltonian was then diagonalized (with the con- 
fining potential omitted) with this augmented basis yield- 
ing a set of positive energy pseudo-states. The phase 
shifts were derived by a least squares fit to the overlap of 
the target and projectile wave functions with the pseudo- 
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states The overlap function, C(tq) is denned as 

C(r ) = J d 3 n d 3 r 2 ^ H2 (ri,r 2 )*(r ,ri,r 2 ) . (6) 

The overlap function depends on the distance from the 
inter-nuclear midpoint and the angle, 8q from the inter- 
nuclear axis. Least squares fits to roC(ro) over the finite 
interval, tq G [Ri , i? 2 ] , at fixed values of 9q were made 
to the asymptotic form _Bsin(fcro + Sq). The radial lim- 
its for the fit were chosen as R\ = 18 ao and R2 = 30 
do- This procedure is reminiscent of an earlier method to 
determine molecular phase shifts using discrete functions 
[22! . The lowest energy pseudo-state was at k « 0.006 
ciq . The scattering length was extracted from the phase 
shift using A scat rj — tan(<5)/fc while Z c ff is determined 
from the normalization constant. Table [J gives the scat- 
tering length and Z c g for the lowest energy pseudo-state 
extracted for projections parallel and perpendicular to 
the inter-nuclear axis. 

An alternate estimate of the scattering length was 
made from the energy. The evenly tempered asymptotic 
positron basis was diagonalized for a zero potential. This 
basis can be regarded as defining a soft-sided box 
The effective radius of this box can be estimated from 
the lowest energy V = state, and the radius allows the 
scattering length and Z e s to be determined (l9j |. These 
are designated in Table Q] as Ascat.SB and Z e ff,SB- The 
methods used to estimate the scattering length do not 
take long range polarization and quadrupole interactions 
into account past ro ~ 24 ao. Subsidiary calculations 
suggest an underestimation of |A scat | by about 5%. 
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FIG. 1: The scattering length (in ao) as a function of inter- 
nuclear distance, R for positron scattering from H2. 

The scattering length in Table U becomes increasingly 
negative as the dimension of the basis increased. This 
is expected on physical grounds. Comparison between 
^-scat, II ) Acat,j_, and Ajcat^B and the Z e ff,||,j_,SB values 
reveals the extent to which the low energy scattering pa- 
rameters are largely unaffected by the aspherical poten- 
tial. The overall variations between the values of Z c g and 



^4scat are about 1%. The calculations at this energy are 
equivalent to the H2 molecule being its lowest rovibra- 
tional level. It must be kept in mind that our calculation 
is for a fixed axially-symmetric target, while a non-Born- 
Oppenheimer calculation would treat the H2 system as a 
spherically symmetric system. 

The UN group had previously used the method of mod- 
els within the Kohn variational method to determine the 
low energy Z e g. The value listed in Table 1 is taken from 
the calculations labelled "ii" in Table 4 of 0. This gave 
a Z e g of 10.3. A Kohn variational calculation which ex- 
plicitly included the H2 wave function was very recently 
reported by the UN group ll|. The result given in Ta- 
ble [T] used a H2 wave function which gave 99.7% of the 

f2 B) 

correlation energy and were taken from the ' curves 
in Figures 7 and 8 of [ll|. Some UN method of models 
calculations published while the present letter was under 
review gave Z e s = 13.5 [Hj]. The same article also gave 
a Z e g w 10 with an explicit H2 wave function and the 
UN group did not make a clear statement about which 
result should be preferred [l2j]. 
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FIG. 2: The close to zero energy, Z e ff, as a function of inter- 
nuclear distance, R. The cross indicates the location of the 
R = 1.448 ao experimental values listed in Table HI 

Calculations have also been performed at a series of 
inter- nuclear separations between 1.0 ao and 4.4 ao. The 
scattering length as a function of inter-nuclear separation 
is shown in Figurc[T]whilc the zero-energy Z e f{ is depicted 
in Figure O Z e g for the vibrational ground state was 
estimated by assuming the linear form Z e g(R) w Z + 
Z\R. The Z c ff for the vibrational ground state is then 
computed by evaluating Z c s at the mean inter-nuclear 
distance, (R). Computing Z e g at (R) = 1.448 ao [23[ 
gives (Zoff)vib = 15.72. The scattering length for the 
vibrational ground state was estimated at —2.74 an. 

Experimental Z eS values of 14.7(2) [|, 16.02(8) Q and 
14.61(14) [| have been measured. The differences appear 
to be related to variations in Z e s with gas density for 
reasons that are not known Q . The present calculation is 
compatible with experiment when consideration is given 
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to the uncertainties in the experimental analysis. The 
traditional gas phase positron annihilation experiments 
simply inject high energy positrons into the gas and rely 
on the assumptions that the positrons are thermalized 
and no other processes are occurring when the lifetime 
spectrum is measured. 

The zero energy vibrational Z c g still needs to be con- 
verted to thermal energies. A rough estimate of the size 
of the correction can be made by using an approximate 
form for the energy dependence of Z e g [24|, e.g. 

7 ,,s Zcff(O) 

Zcs{k) = i + (A scatk y ' (7) 

Application of this result with a scattering length of —2.7 
ao suggests a 1.5% reduction in the annihilation param- 
eter at thermal energies to a value of 15.5. 

The scattering length implies a zero energy cross sec- 
tion of g(0) w 30 nay. A recent experiment by the Trcnto 
group [25| had a cross section of 8.3 ita^ at k « 0.086 a^ 1 . 
The experimental cross section is absolutely incompati- 
ble with the present scattering length and that of the 
UN group @. Improving the quality of the CVM wave 
function would only lead to the magnitude of the scatter- 
ing length increasing, thus lea ding to larger discrepancies 
with the Trento cross section |25| . 

The scattering length shows a tendency to increase in 
magnitude as the inter-nuclear separation is increased 
and a virtual state is formed around R w 3.4 ao. The 
maximum scattering length is —13.0 ao at R = 3.4 ao- 



The peaking of Z e g around 3.4 ao is expected since it 
is known that a large scattering length leads to a large 
threshold Z e g [24[ . The large scattering length was a sur- 
prise. However it is known that the critical value for an 
electric quadrupole to bind a charged particle is 2.4 ea\ 
[2(|. The quadrupole moment of H2 increases from 0.91 
eao at R — 1.4 ao before reaching a maximum value of 
2.03 eao at i? = 3-0 a [27]. We speculate that the large 
increase in scattering length can be understood in terms 
of the larger quadrupole moment. The recent method of 
models calculation by the UN group exhibited a qualita- 
tively similar variation of Z e g versus R [l2T ] . 

While the present calculation was performed un- 
der the fixed nucleus approximation, it represents the 
first description with an unrestricted treatment of the 
positron/electron interactions in the e + -H2 collision sys- 
tem. The strong increase in Z c g and A scat with increasing 
inter-nuclear distance due to virtual state formation at 
R w 3.4 a was totally unexpected. One of the most sig- 
nificant methodological aspects was the ease with which 
the inner region wave function was generated. Using the 
present e + -H2 wave function within a more formal scat- 
tering framework, such as the Kohn variational method, 
would require substantial development work, but this 
would involve the application of known procedures and 
would be straightforward. 

This work was supported under the Australian Re- 
search Council's Discovery Program (project number 
0665020). 
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